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Summary: 

We study the influence of disorder and randomly distributed impurities 
on the properties of correlated antiferromagnets. To this end the Hubbard 
model with (i) random potentials, (ii) random hopping elements, and (iii) 
randomly distributed values of interaction is treated using quantum Monte 
Carlo and dynamical mean-field theory. In cases (i) and (iii) weak disorder 
can lead to an enhancement of antiferromagnetic (AF) order: in case (i) 
by a disorder-induced derealization, in case (iii) by binding of free carri- 
ers at the impurities. For strong disorder or large impurity concentration 
antiferromagnetism is eventually destroyed. Random hopping leaves the lo- 
cal moment stable but AF order is suppressed by local singlet formation. 
Random potentials induce impurity states within the charge gap until it 
eventually closes. Impurities with weak interaction values shift the Hub- 
bard gap to a density off half-filling. In both cases an antiferromagnetic 
phase without charge gap is observed. 



1 Introduction 



Antiferromagnetic spin correlations are present in many strongly correlated elec- 
tron systems, notably the prototype Mott insulators NiO and V2O3, the parent 
compounds of HTSC cuprates, and heavy fermion systems such as YbP, U2Zni7, 
and many others. Many of those systems are intrinsically disordered, in partic- 
ular upon additional homo- or heterovalent doping. The influence of impurity 
doping on antiferromagnetic (AF) order and electronic properties has recently 
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been studied in a variety of systems. Doping with static scatterers like nonmag- 
netic impurities usually weakens antifcrromagnetic order, a prominent example 
being Zn doping in YBa2Cu306 [fj]. In spin chains (CuGeOs) 0] and ladder com- 
pounds (SrCu20s) |}| on the other hand, doping with magnetic and nonmagnetic 
impurities can induce AF order while the pure systems show spin gap behavior. 
Very effective in destroying AF order are mobile carriers, e.g. hole doping in 
Lai_ x Sr x Cu04 Q. The stability of AF order strongly depends on the positions 
of the dopant level. While in the nickel oxides Lai_ x Sr x Ni04 |^] and Nii_ x Li x O 
H holes are supposed to be localized, in the cuprate Lai_ x Sr x CuC"4 the hole 
level lies in the valence band leading to mobile scatterers. As a result, AF order 
is stable in Lai_ x Sr x Ni04 up to x = 0.5, but is destroyed in the cuprate already 
at approx. 5% Sr doping. 

In the present paper we will study the influence of disorder on AF order and the 
Mott band gap in correlated antifcrromagnets. We employ the Hubbard model 
in the presence of different types of disorder. While the (disordered) Hubbard 
model is certainly far too simple to describe real materials it already contains 
very rich physics including local moment formation, magnetic ordering, Mott- 
Hubbard transition, and Anderson localization. On the other hand, the interplay 
of disorder and interactions in electronic systems belongs to the most difficult 
problems in physics, and reliable results within simple models are still very 
desirable. The problem has been investigated in the past by a variety of methods, 
including field theoretical approaches M, renormalization group treatments |§, 
[|, unrestricted Hartree-Fock Jl(], [llj , dynamical mean-field theory (DMFT) [ 12 , 
L| ||, quantum Monte Carlo (QMC) @, H, 0, and several more (see [if] 
for a review). Here we give an overview of results obtained by QMC and DMFT 
concentrating on the AF phase diagram and the Mott gap. 

We consider the following Hubbard Hamiltonian: 

H = ^(e, - ^)n ta + M4tCj<x + h.c.) + Y Ui(hif - i)(n a - i). (1.1) 

ia (*i) <T ' 

In principle all parameters €i,tij,Ui can be randomly distributed. The precise 
definition of the different disorder types studied in this paper will be given in the 
following sections. The average t = (tij) sets our energy scale. We will restrict the 
hopping tij to nearest-neighbors hence not allowing for frustration. Longer range, 
random hopping amplitudes will be important in the modeling of amorphous 
materials such as doped semiconductors [|l8[ and are not considered in the 
present work. 



2 Methods 
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2.1 Determinant quantum Monte Carlo (d=2) 

We use a finite temperature determinant quantum Monte Carlo method ppf to 
obtain approximation-free results for finite lattices. The algorithm is based on 
a mapping of the interacting electron problem onto a d + 1 dimensional quasi- 
classical problem using auxiliary Ising-type spins. It provides for calculating 
thermal averages of observables, A, at a temperature T = 1/(3, 

(A) = ^. 2.2) 

The phase space sampling over the auxiliary field configurations is performed 
using Monte Carlo techniques. The weight of a configuration is proportional to 
a product of two determinants, one for each electron spin species. In the case of 
half-filling without random potentials, i.e. e,- — fx = 0, on a bipartite lattice the 
determinants always have the same sign, hence their product is always positive 
semi-definite, which can be shown by particle-hole transformation of one spin 
species [cij. — > (— l^c-i]. In the general situation that the product can become 
negative the algorithm can still be employed in principle. However, the signal 
to noise ratio decreases exponentially with systems size, inverse temperature, 
and interaction, putting severe restrictions to the applicability of the method. 
This so-called "minus-sign problem" is a general obstacle for all exact fermionic 
Monte Carlo methods as well as for spin-systems in the presence of frustration. 
Even without the minus-sign problem the computational effort is large because 
the computer time grows cubically with system size N, restricting N to the order 
of 100 on present supercomputers. 

In the case of disorder all observables have to be averaged over the (frozen) 
disorder configurations. Because of the computational effort we restrict ourselves 
to two dimensional lattices with linear size up to L x = 10 which often allows a 
reliable finite size scaling. Since we are interested in AF ordering we calculate the 
magnetic correlation functions C(l) and their Fourier transforms, the magnetic 
structure factors 5(q), 

C(l) = l^(m jfflj+1 ), 5(q) = ^C(l)e^. (2.3) 
j l 

In particular the AF structure factor S(tt,it) is used to obtain the ground 
state sublattice magnetization M by a finite size scaling Ansatz according to 
spinwave theory 
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2.2 Dynamical mean-field theory (limit of infinite dimensions) 

The dynamical mean-field theory |p2| , |23| is a local approximation in which the 
self energy becomes site diagonal, or momentum independent: 

E y (u;) = <%£(w), E(k,w) = S(o>). (2.5) 

The one-particle Green function G(k, cj) can hence be obtained from the non- 
interaction Green function G°(k, u>) by G(k, u>) — G°(k,uj — S(w)), and the 
local Green function is given by Gu(uj) = I/A^ k G(k, w) This does not imply 
a simple shift of energies, like in traditional mean-field theories (e.g. Hartree- 
Fock), because S remains dynamical, i.e. frequency dependent, preserving local 
quantum fluctuations. The local approximation becomes exact in the limit of 
infinite spatial dimensionality and maps the interacting lattice model onto a 
self-consistent single impurity model like, for example, the Wolff model p4| ]: 

H Wo m = ^ £k "- k + u ( 7 " l oi - gX^oj. - ^) + e ("-ot + noi). (2.6) 

kcr 

Here the one-particle energies have to be defined such that the non-inter- 
acting local Green function of the Wolff model fulfills G^ olff = (Gu — E) ~~ 1 . In the 
self-consistent solution the local (interacting) Green function of the Wolff model 
has to be equal to Gu. In the presence of disorder one has to average over all 
possible values of U or e, respectively. This type of local averaging is equivalent 
to the "coherent potential approximation" , well known from investigations of 
disordered alloys. 

While the self-consistency is rather easily reached by iteration the solution of 
the single impurity problem is the hard part. There exists no analytic solution 
and different numerical and approximative techniques have been employed (see, 
e.g., Q). Here we again use auxiliary field QMC [^5), an algorithm quite similar 
to the one for finite dimensional lattices sketched above. The computer time 
grows like L 3 where the number of Matsubara frequencies, L oc (3. Fortunately, 
QMC for the single band model is free from the minus-sign problem. 

In the following the non-interacting DOS is chosen as a semi-elliptic model 
DOS with bandwidth=8, equal to the d — 2 tight binding bandwidth for t = 1. 
A typical quantity under consideration is the staggered magnetic susceptibil- 
ity xaf whose divergence signals the transition to an AF ordered state. One 
can also extend the DMFT equations to the ordered phase to obtain spin and 
sublattice dependent electron densities and the sublattice magnetization M. The 
one-particle density of states (DOS) is obtained by analytical continuation of the 
imaginary time Green function using the Maximum Entropy method. For details 
of the algorithm, the implementation of disorder averages and determination of 
expectation values sec |23|, Q . 
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3 Random potentials 
3.1 Local moment quenching 

Random potentials are the most frequently studied type of disorder in the con- 
text of Anderson localization. Contrary to the Hubbard interaction which at 
half-filling favors single occupation on each site, different local potentials lead 
to different occupations and hence to a quenching of local magnetic moments 
on sites with large absolute value of the local potential. This is seen in Fig. la 
where the average local moment squared, m 2 , is plotted versus disorder strength 
for a flat distribution of e values with width A. It is also shown that the spin- 
spin-correlations go in parallel with m 2 . For a large width of random potentials 
this type of disorder is apparently very effective in the destruction of magnetic 
order. To study the behavior of the charge gap, the electronic compressibility, 
K = dn/d/j,, is calculated as a function of A (Fig. lb). While disorder decreases 
k in the non-interacting case, k is enhanced by disorder at finite U . The reason 
is the introduction of states within the AF charge gap as will be discussed be- 
low. As mentioned above, random potentials break the particle-hole symmetry 
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Figure 1 (a) Local magnetic moment C(0, 0), spin-spin correlation function at 
largest possible distance (2,2), and AF structure factor S(ir,iv) as a function of 
disorder strength for a constant distribution of random potentials (in all cases 
the values at U = A = are subtracted), (b) compressibility k vs. A; disorder 
reduces k in the non- interacting case and enhances it at U = 4 |16|. 
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leading to a minus-sign problem even at half filling. This is the reason why only 
small lattices (4 x 4 in Fig. 1) can be studied for this type of disorder. 



3.2 Disorder-enhanced derealization 



Consider the situation deep in the AF phase, i.e. with a large staggered moment. 
A high value of the local potential on a given site reduces the potential barrier 
for a majority spin electron to tunnel to a neighboring sites. The electron is 
hence delocalized and the local magnetic moment is reduced. A low potential on 
the other hand cannot significantly further localize the majority spin because it 
is already almost saturated. This asymmetry of the localizing and delocalizing 
effects of random potentials is depicted in the DMFT results, Fig. 2. The ma- 
jority spin (|) is strongly reduced for e > but almost unchanged for e < 0. 
Note that the total local density monotonically decreases with e, and that the 
net magnetization (nj — n±) is reduced for any |e| > 0. 



CD 




Figure 2 Local spjn resolved densities in DMFT as a function of the random 
potential value a JM. Width of disorder distribution is A = 8(=bandwidth). 

The fact that the derealization is not compensated can be expressed in terms 
of an enhanced effective hopping parameter t c g. In the case (U — A) ;§> t the 
effective i e ff can be estimated by the mapping onto the AF Heisenberg model: 
the AF exchange arises from virtual hopping of an electron with spin a to a 
neighboring site occupied with an electron of spin — a. The exchange energy for 
this and the opposite process is 
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Figure 3 AF (T — (/) phase diagram in DMFT for a bimodal random potential 
distribution. Energies in units of the half bandwidth. Jl4] 



Jij is hence always enhanced from the pure case (Jo = 2t 2 /U), and the average 
exchange becomes J av = Jq[1 + A(A/C/) 2 ] where A depends only on the shape 
of the disorder distribution (e.g., A = 1/12 for the constant, A = 1/2 for a 
binary distribution). Furthermore, one can show that the form of the magnon 
propagator remains unchanged in second order in A/U p6| , but magnons are 
stiffened by a factor of (1 + X(A/U) 2 ). The effective hopping in strong coupling 
can thus be expressed as t e g = t[l + A(A/J7) 2 ] 1 / 2 . In dimensions d > 2 one 
expects the Neel temperature Tjv to be proportional to J and hence Tjv, too, 
should be enhanced by weak disorder. 

While disorder-enhanced derealization stabilizes AF order at strong coupling 
it suppresses it at weak coupling. For small U the AF state is rather "spin 
density wave like" and an enhanced kinetic energy tends to weaken AF order 
and to reduce TV- In addition, random potentials destroy the perfect nesting 
instability responsible for AF order at small U. 

The T vs. U phase diagram within DMFT (Fig. 3) summarizes and confirms 
the above considerations: For small U, Tn is reduced and eventually vanishes 
when A becomes roughly equal to U. At large U, however, the Tjv curves for 
different A cross each other, i.e. at a given value of U, Tn increases with dis- 
order. The opposite effects of disorder on Tjv depending on U are due to the 
non-monotonic behavior of the function Tjv(/7). If the slope oiTfj(U) is positive 
(negative) an effectively reduced ratio U/t e g leads to a suppression (enhance- 
ment) of Tpj. 
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3.3 Closing of the charge gap 

Local potentials can induce electronic states in the charge gap. As also found 
experimentally the positions of the impurity states are crucial for the stability 
of AF order with respect to carrier doping. Disorder-induced gap states will 
reduce the size of the charge gap. We now want to determine the critical dis- 
order strength at which the gap vanishes. First we treat the disorder within 
the T- matrix approximation which becomes exact for a single impurity |2^, [2^] . 
Since the local host Green function of the correlated, not-disordered problem 
is unknown we approximate it by the Green function of the AF Hartree-Fock 
solution: 

[g% t = (l/N) 5> - aD)/{^ El). (3.8) 

k 

Here 2D = mU is the Hubbard energy gap in the pure AF, = yj D 2 + e£ 
is the AF band energy, and D is obtained from the self-consistency condition 
(l/N) X)k(^^ k ) _1 = U^ 1 . With this local host Green function we can calculate 
the location of impurity-induced states from the poles in the T- matrix, T a (ui) = 
Ci/(1 — Ci{g°(uj)]ii) . We consider a constant distribution of random potentials 
between ±A/2. There are disorder-induced states within the gap if |[<?°]m| > 
2/ A. The energy D up to which states are formed within the gap is given by 
[g®(— D)\n = 2/ A. The remaining charge gap 2D is plotted versus A in Fig. 4. 
The decrease is almost linear and 2D vanishes close to A = U. Also shown is 
the result from a numerical unrestricted Hartree-Fock (UHF) analysis. In this 
approach, the HF Hamiltonian on a finite lattice is numerically (self-consistently) 
diagonalized, so that disorder is treated exactly O, |26|. The energy gap is 
obtained from the energy difference between the lowest energy state of the upper 
Hubbard band and the highest energy state of the lower Hubbard band. Averages 
are taking over 100 disorder realizations on a 10 x 10 lattice. The saturation of 
the gap at A/U ~ 1 is due to the finite system size. The agreement with the T- 
matrix approach is excellent for the present interaction value U = lOt. Deviations 
from the T-matrix approach are more pronounced at lower interaction strengths 
where the fermion states are more extended. 

The closing of the charge gap is also observed within the DMFT approach. 
Fig. 5a shows the density of states for several disorder values. Note that all 
spectra shown are within the AF ordered phase, i.e. AF order is much more 
stable than the charge gap which closes about A = U. A linear reduction of the 
charge gap is also observed in the case of a binary disorder distribution (Fig. 5b). 
For this stronger type of disorder both the charge gap and the AF order vanish 
near A = U. 
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(a) (b) 

Figure 5 (a) Density of states (DOS) for a constant disorder distribution of 
different width in DMFT |26|; all spectra are within the AF phase, (b) Charge 
gap vs. A in DMFT for theDimodal distribution of random potentials. 



3.4 Spin vacancies 

As discussed in the previous section the one-particle excitation gap decreases 
with A and vanishes near A = U. Upon further increase of A the two Hubbard 
bands will overlap whereby electrons from the highest levels of the lower band 
(with £j > U/2) will be transfered to the lowest levels of the upper band (with 
Ei < —U/2). Those sites will become doubly occupied, their local moment will 



10 



Ulmke et al. 



vanish as observed in Sec. 3.1. In the limit t <C U, A the Hubbard bands at ±17/2 
are nearly flat with width A. The overlap region is ~ (A — U) and the fraction 
of nonmagnetic sites can be estimated by x ~ (A — U)/A. 

In the following this situation is modeled by spin vacancies of concentration 
x which will lead to strong magnon scattering. Such diluted antiferromagnets 
have been widely studied recently mostly within spin diluted Heisenberg models 
pgl |29f and also Hubbard models ^7j. A perturbative analysis in the strong 
coupling limit {0(t 2 /U)) (2?], ||(| yields a softening of magnons by a factor of 
(1 — 2x) in d = 2. Extrapolating to large fractions x, the magnon energy scale 
J = J(l — 2x) hence vanishes only close to the percolation threshold x pcrc ks 0.4. 
In d > 2 the Neel temperature is assumed to be proportional to J and a linear 
decrease of Tjv is thus expected. In experiments, e.g. on Li doped NiO such 
a linear dependence on x is indeed observed with a prefactor of 2.2. In the 
limit d —* oo the percolation threshold is 1.0; hence we can expect AF order 
to persist for arbitrarily large dilution. Within DMFT, too, a linear decrease 
of T/v oc (1 — a;) is observed for the diluted model at large U (Fig. 6a). At 
small U, however, AF order is remarkably robust, Tjy being constant up to 
dilutions close to x = 1 where Tjv eventually drops to zero. This behavior of T/v 
can be explained by a strong enhancement of the local moment density on the 
remaining sites at small U (Fig. 6b). The reason for this enhancement is that 
with reduced average number of nearest neighbors the kinetic energy decreases, 
leading to a stronger localization. At large U the local moments are already 
almost saturated at x = and just cannot be further enhanced. The situation 




Figure 6 (a) Tn as a function of site dilution in DMFT for different values of 
U. Dashed line Hartree-Fock approximation for U = 0.65 (energies in units of the 
half bandwidth), (b) Enhancement of local moments on the remaining fraction 
(1 — x) of the lattice at inverse temperature /3 = 16. H] 
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in the case of vacancies is therefore quite different from the effect of weak site 
disorder where disorder-enhanced delocalization is observed. 



4 Random hopping 

The case of spin vacancies discussed in the previous section can of course be 
regarded as a specific type of randomness in the hopping elements. The more 
generic case of a continuous (flat) distribution of £y G [1 — A/2, 1 + A/2] was 
studied in d = 2 using QMC. Since the hopping is still restricted between nearest 
neighbors on a square lattice particle-hole symmetry is preserved and no minus- 
sign problem occurs at n = 1. Random tij hardly affect the density of local 
moments (Fig. 7a), the slight decrease may be due to the enhanced kinetic 

energy which is oc ./ '< tfj >. Nevertheless, longer-range spin-spin correlations 

are strongly suppressed if A is of the order of t (see Fig. 7a). The finite size 



scaling according to (2.4) yields the AF order parameter (staggered moment) M 
versus A (Fig. 7b). M vanishes at a critical disorder strength of A c w 1.4. We 
propose |l6[ that the phase boundary is determined basically by the variance of 
the AF exchange coupling v — (< Jfj > — < Jy > 2 )/ < Jfj >. AF order persists 
for v < v c w 0.4. This criterion is consistent with the phase boundary of the 
bond-disordered AF Heisenberg model with a bimodal distribution of Jy |3lj]. 
The reason for the vanishing of AF order for this type of disorder is supposedly 
the formation of local singlets. Such singlets will form first on the strongest 
bonds and will leave some spins which are weakly coupled to their neighbors 
unpaired. Those "free" spins are expected to give a Curie-like contribution to 
the susceptibility as observed in doped semiconductors Jll| fl9). The numerical 
results indeed show a strong enhancement of the uniform susceptibility in the 
disordered case 1l6f . 



5 Impurities with weak local interactions 

The Hubbard model at n — 1 exhibits both AF order and the Mott-Hubbard 
metal- insulator transition (MIT) . It is important to note that both effects are in 
principle independent. The Mott-Hubbard MIT occurs at intermediate interac- 
tion (U ~ bandwidth) while AF order can set in at arbitrarily small U due to the 
Fermi surface nesting instability. It is argued that the MIT at T — represents 
a quantum critical point which is however concealed by the low temperature AF 
phase |32Tl , and attempts are made to suppress the AF phase to very low or zero 
T by different types of frustration A different possibility to separate AF 
order and the MIT is to shift the filling at which the Mott-Hubbard MIT occurs 
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Figure 7 (a) Local moment and spin-spin correlations in d — 2 for random 
hopping (same quantities as in Fig. 1). (b) Staggered magnetization M vs. A as 
obtained by finite size scaling Jl^ ]. 



away from n — 1 by introducing impurities with a low (or zero) local interaction 
with a concentration /. Such impurity sites can be doubly occupied without the 
cost of the local repulsion and hence the MIT is expected at a density n = 1 + / . 



5.1 Mott-Hubbard gap 

The Hubbard model with a bimodal distribution of U values (U = on a fraction 
/ of the lattice, and U = 8t on the remaining sites) was investigated by QMC on 
square lattices and in DMFT. Note that for technical reasons (minus-sign prob- 



lem) the choice of the [/—disorder in model (1.1) preserves particle-hole sym- 
metry. This corresponds to different chemical potentials on the two constituents 
such that at n = 1 the electronic density is homogeneous, i.e. independent of 
Ui. To detect charge gaps the electron density n is plotted versus the chemical 
potential fi (Fig. 8a). As expected the gap moves to a density off half filling close 
to 1 + /. This agrees with the kinetic energy which shows minima at the cor- 
responding densities (Fig. 8b). At a closer look one observes an additional gap 
at ri = 1 which is due to the doubling of the unit cell in the AF ordered state. 
The gaps at ti = 1 and n = 1 + / can also be detected by the one-particle spec- 
trum, obtained in DMFT (see 0|). Between densities 1 and 1 + / the systems 
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(a) (b) 

Figure 8 (a) Electron density n vs. chemical potential n for different values 
of U = impurity concentration /. n(/i) has plateaux at n = 1 and close to 
n = 1 + /. (b) Kinetic energy (k) vs. average density for f — (open triangles) 
and / = 0.111 (solid triangles), (k) shows minima at n = 1 + / Em. 



is apparently compressible. 

5.2 Antiferromagnetic order 

Figure 9a shows the staggered moment M extrapolated to the thermodynamic 
limit in d = 2 at n = 1 as a function of impurity concentration /. For small / 
AF order is robust leading to the charge gap at n = 1. AF order vanishes at 
f c « 0.45, i.e. close to the percolation threshold / per c ~ 0.4. 

Off half-filling the finite size scaling is no longer possible due to the minus-sign 
problem and the AF phase boundary is obtained within DMFT only (shown 
in Fig. 9b at T = 1/8). It is found that the U = impurities can induce 
AF order at densities for which the clean model is disordered. In the clean 
model (/ = 0) additional electrons are free to move and hence very effective in 
destroying long range order. U — sites provide localizing centers which are 
energetically favorable for the additional carriers. This is why the AF phase 
extends to larger dopings at finite (but small) values of /. Since the mechanism 
of localizing mobile dopants is observed in d = 2, too, we expect the enlargement 
of the AF phase to be present in d — 2 at T — as well, in spite of the fact that 
in the clean model in d = 2 the critical doping is supposed to be zero. Eventually 
at larger fractions / the AF phase shrinks, and the critical density approaches 



14 



Ulmke et al. 



U = 8t T=1/8, U=8 




(a) (b) 

Figure 9 (a) Staggered magnetization in d = 2 as a function of impurity 
concentration /. (b) f — n phase diagram at temperature T — 1/8 within DMFT. 
For small / the AF phase is stabilized against doping. |0| 



1.0 for f c « 0.6 (/ c « 0.75 in the ground state jnj). For even larger values of / 
there is no AF order even at n = 1. 

Both the stabilization of AF order and the shift of the Mott gap to higher 
densities result from the localization of carriers at the {7 = impurities. AF 
order is only destroyed when the density on the U = sites saturates, i.e. it is 
stable across the Mott-Hubbard MIT. Hence the separation of Mott-Hubbard 
MIT and AF order is not present in the U = impurity model, at least not for 
the present parameter values. 



6 Summary and Conclusion 

In this paper we discussed different types of disorder in correlated antiferromag- 
nets and presented results obtained mostly by quantum Monte Carlo simulations 
in d — 2 and within dynamical mean- field theory (d = oc). 

Different mechanisms were identified by which disorder can enhance antiferro- 
magnetic order: i) disorder-enhanced derealization at strong coupling in the case 
of weak disorder in the chemical potentials and ii) localization of surplus carriers 
in the case of impurities with weak local interaction. In both cases compressible 
antifcrromagnetic phases are observed. To determine if the gapless AF phase 
is metallic requires the calculation of transport properties which is presently in 
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progress. 

Quantum Monte Carlo simulations of electronic tight binding models are just 
making the transition from addressing rather abstract issues of correlation effects 
to making contact with real experiments. One important feature in this respect is 
the treatment of intrinsic disorder. For a quantitative description of experiments, 
however, the inclusion of a realistic bandstructure is mandatory. Here the DMFT 
will be particularly helpful because it allows for the treatment of multiband 
models in a wider range of parameters and in the thermodynamic limit. 
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